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We study geodetic lines on a surface generated by a small deformation of the standard 2D-spliere. 
We construct an auxiliary hamiltonian system with the view of describing geodetic coils and almost 
closed geodesies, by using the fact that loops of the coil can be well approximated by great circles 
of the sphere. The phase space of the auxiliary system is determined by the graph generated by 
separatrixes of its solutions, the vertices of the graph corresponding to almost closed geodesies 
and the edges to the geodetic coils joining them. Topological types of the graph depend on the 
parameters determining the deformation. Using the method of averaging in conjunction with the 
computer modelling of the auxiliary system, we obtain a fairly detailed visualization of geodesies on 
the deformed sphere. 
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I. GEODETIC COILS 
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Geodesic lines on a surface can be considered either as straight lines with respect to a Riemann metric, or trajectories 
of a particle of mass m moving freely on the surface. The second approach allows for using the methods of analytical 
mechanics, and has drawn considerable attention, flj. It should be noted that the general solution of the problem 
of geodetic lines is known only for certain special cases, e.g. ellipsoid, 0,01 and the topology, or analysis situs, of 




FIG. 1: Two coils I, II of a geodesic on ellipsoid; vectors Ni and Nu are the normals to the planes of great circles approximating 
the coils 
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FIG. 2: Ellipsoid with main axes 1.01, 1.02, 1.03; the geodesic, corresponding to a separatrixe of the auxiliary system, leaves a 
saddle point; the initial point (-0.0117, 0.0001, -1.0299); the initial velocity (-1.9416, 0.0194, 0.0228.) 




FIG. 3: Ellipsoid with main axes 1.01, 1.02, 1.03; the geodesic, corresponding to a separatrixe of the auxiliary system, arrives 
at a saddle point; the initial point (0.7416, 0.6169, 0.3176); the initial velocity (0.4728, -1.2782, 1.3832) 



geodesies on a surface needs specific studying. In this paper we use asymptotic methods and in particular focus our 
attention on geodetic lines that are closed to within the first order of perturbation theory. 

The central idea relies on the circumstance that in case of an ellipsoid not differing substantially from a sphere, 
the great circles of the latter may serve a good approximation to the ellipsoid's geodesies, if they are short enough, 
and one can visualize the geodesies as winding up in coils; loops or rings of the coil corresponding to great circles of 
the sphere, see FIG^ Hence approximating the successive rings by great circles, we may describe the change in the 
position of the rings by the motion of a great circle, see FIG^ which in its turn is determined by the normal vector L 
of the plane cutting the sphere along the great circle. To cast this picture in a more quantitative form we may use the 
fact that the normal vector L is the angular momentum of the particle moving along the geodetic line. This, we obtain 
the advantage of allowing for the use of asymptotic methods, in particular the averaging. Our main instrument is an 
auxiliary hamiltonian system describing the asymptotic motion of L. It can be easily solved, so that its stationary 
solutions correspond to closed planar geodesies within the limit of precision provided by the averaging method , 
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see FIG|21 We can get some insight into the behaviour of geodesies by considering solutions joining the stationary 
solutions, that is separaterixes, which, from a purely geometrical point of view, are a kind of bridges between closed 
geodesies. In the case of two dimensional torus in three dimensional spcae this phenomenon had been observed by 
Yu.S.Volkov, 0]. Generally, a separaterixe geodetic coil does not coincide anywhere with a closed geodesic but comes 
infinitely close to it, and thus resembles a limit cycle, see FIG|2 The set of stationary points and separaterixes of the 
auxiliary system forms a net, or graph, which, as is shown in Section II can be realized on projective plane. Thus, we 
obtain a topological classification of geodetic coils, and find that the number of their topological types is finite. 

These arguments are valid even for surfaces having a more general form than that of the ellipsoid. Studying the 
specific case of the deformation determined by fourth order terms involves the surface being substantially different 
from the ellipsoid as regards its differential geometry, as well as the algebraic structure of its equations, and hence we 
come across considerable difficulties by attempting a direct investigation of the geodesies. The asymptotic approach 
outlined above has the advantage of getting round this problem. 



II. AVERAGED EQUATIONS OF GEODESICS 

The equations determining geodesies on a surface can be cast in the form of the equation, 
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where f{x) is the right-hand side of the constraint f{x) = determining the surface. The Lagrangian multiplier can 
be found explicitly , so that the equation of motion, in the form that does not involve A, reads 



'dipY dx- 
dx 

In this paper we consider surfaces that do not differ substantially from sphere, and in assume that their equations be 
of the form 



i 

where Si are small. Then equations lU read 

The equations given above are more tractable than the usual ones employing the Ghristoffel symbols and an explicit 
parametrization of sphere, so that one may prefer them for the needs of numerical simulation, as is done in this paper. 
With the view of obtaining a qualitative description of geodesies, we shall consider the angular momentum 

L = r X p, 



Its components satisfy the equations, which follow from Q 
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Even though the equations given above are exact, in the sense that they do not involve any approximation and do 
not use the being small, their treatment still need further refining, and this will be done with the help of the 
method of averaging. Generally, the approach relies on studying the evolution equations for integrals of motion of 
the unperturbed system, i.e. in our case the normals to the planes of the large circles, with respect to the basic 
periodic solution of the latter. The averaging serves as a filter separating the main regular part of the solution from 
the oscillating one caused by small terms considered as perturbation, see 0. 

We shall write the basic equation for the particle's motion on the sphere of unit radius in the form 



X — cos{ijjt + 6')ei + cos{ujt + 9)e2 
vectors ei, 6*2, 63 determined by the equations: 
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(0,^3,-^2) 



— L1L2, L1L3) 



6*3 = — (Li,L2,i3)- 



The angular velocity uj is given by the equation uj"^ — x — L^, valid to within the first order of perturbation. Here 
Li,L2,Lj, are coordinates of the normal to the plane of the great circle determining the solution, i.e. the angular 
momentum. 

Let us turn to the exact equations for the angular momentum Q . With the help of the equations given above and 
neglecting terms of the second, and higher, order in the Si, we can transform equations Q in the form 
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{L 



L1L2 



cos{ujt + 9)L3 + sin(a;i + i 



cos(wt + 0){-L2) + sin(wt + 6) 



cos{ujt + 6')(— L2) + sm{ujt + I 
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cos(a;t + d)L3 + sm{ujt - 
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L2 = 



2L^e, 



cos{ujt + 9){—L2) + sin(a;i + 
cos{ujt + 9) -0 + sin(tjt + 6)- 



^2 ~ ^3 
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cos{ujt + 9) - + sin{ujt + 9)- 



-L\ — Li 



cos{Lot + 9)L^ + su\{Lut + 9) 



L1L2 
L 



cos((jjt + 9)L3 + s'm{ujt + 9) 



L1L2 
L 



cos(ajt + 6*) • + sin(a)t + 1 



It should be noted that the right-hand sides of Eqs.Q comprise terms oscillating in time and terms that vary slowly. 
The situation can be treated within the framework of the averaging method, 5], that is on neglecting the oscillatory 
terms we obtain the averaged equations for the angular momentum 



^1 = \ ^[{e,-e2)Ll+e,Ll-e2L'^ 



Ln = - 



3 L^Li 
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[-£3i? + (£l- £3)^2 +£1^3] 



(6) 



4 = \ ^[e2Ll-e,Ll + {e2-e,)L'^ 



It is worth noting that Eq. Q have the Hamiltonian form determined by the usual Poisson brackets for the angular 
momentum, 0, Q, 

{Li, Lj} = ^ SijkLk, 

k 

and the Hamiltonian 
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This circumstance is particularly interesting because, usually, the averaging procedure is not compatible with Hamil- 
tonian structure. 

We may infer from the two integrals of motion, L^ and H, that they admit of an explicit exact solution that can 
be cast in the form of the equation 



iL2 



dLi 



in which 
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-ei(e2 + £3)iL' - LfY - S2S3iL' + Lf)' + —{e2 + e3)L'H. 



The exact solution is not particularly practical. In this paper we shall not use it, but employ topological consider- 
ations and numerical simulation so as to understand the dynamics of L and, consequently, the behaviour of geodetic 
coils. The important point is considering the stationary solutions to Eqs.lQ for which the right-hand sides turn out 
to be zero, and split into three parts SI, S2 and S3, determined by conditions on as follows. 



SI No algebraic constraints imposed on sf. 
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a. LiQ = 0, L20 = 0, L30 ^ 0; 

b. Lio = 0, L20 ^ 0, L30 = 0; 

c. Lio 7^ 0, L20 = 0, L30 = 0; 

S2 The constraints on L relaxed and linear constraints imposed on sf. 



a. 


Lio = 0, 


^20 




^30 0, 


£3^20 


- £2^30 
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b. 


L20 = 0, 


-^30 




^10 ^ 0, 


£1-^30 


- £3-^10 


= 


c. 


L30 - 0, 


Lio 


^0, 


^20 7^ 0, 


£2^10 


- £1-^20 


= 



S3 Vector L subject to Lio 7^ 0, L20 7^ 0, L30 7^ and the quadratic constraints imposed on ef. 

t2 t2 t2 
^10 _ ^10 _ ^10 

£1 £2 — £2 £3 + £3 £1 £1 £2 + £2 £3 ~ £3 £1 ~£i £2 + £2 £3 + £3 £1 

It is worth noting that equations S2 involve the fulfilment of the inequalities £2£3 > 0, £361 > 0, and £122 > for 
cases S2.a, S2.b, S2.c, respectively, whereas equations S3 involve 

£i£2-£2£3+£3£i > 

£l£2+£2£3-£3£l > 
-£l£2 +£2£3 +£3£l > 

Linearizing Eqs.(|HJ at the stationary solutions and, considering small fluctuations of L round them, we may study 
their stability, which turns out to be determined by the requirements 

51 a. £i£2 > 0; 
b. £2£3 > 0; 
C. £3£i > 0. 

52 a. £i£2 - £2£3 + £3£i < 0; 
b. £i£2 + £2£3 - £3£i < 0; 
C. -£i£2 + £2£3 + £3£l < 0. 

53 any Si . 

We may put these equations in a more graphic form using the integral L^ = const, and consider the motion of L on 
a sphere of fixed radius, the integral of energy H taking appropriate values. Then the stable solutions are fixed points 
as regards Eqs.@, the stable and the unstable points are foci and saddle points, respectfully, the separaterixes being 
lines joining the fixed points. Together, they generate a graph on the sphere, having the fixed points as vertices and 
the separaterixes as edges. It is important that the separaterixes, i.e. the edges of the graph, are oriented according 
to the time, so that the graph is the oriented one, and invariant with respect to the symmetry R — > — i?, t —t. 

The pictures of the separatrixe net on the sphere are rather difficult for visualizing, and therefore there is a need for 
some means which are easier to implement. We shall employ the familiar construction of the projective plane, which 
runs as follows. Given a pair of twin points R and —R belonging to the sphere we shall choose the representative of 
the pair determined by the condition R3 > 0; this enables us to consider only the upper part of the sphere. Next 
take the projection of the upper hemi-sphere on x-y plane along axes-z. Thus, a pair of antipodal twins obtains a 
representative in the disk, i.e. a point inside the disk or a pair of antipodal points at the boundary. The separatrize 
nets constructed in this way on the projective plain, are shown in FIGslQ] — I?)). It is important that the equations of 
motion are invariant under the transformations Pi, i = 1, 2, 3 

P^■■ t-. -t, L, -L, (9) 

which generate automorphisms of the graph. 
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FIG. 4: Separatrix net corresponding to the graph of Type I on the sphere; the twin points correspond to the symmetry given 
by Eq.© 




FIG. 5: Diagramme of the projective plain by means of a disk with identified antipodal points at the boundary circle. 



From the fact that the transformation 

P: t^-t, L,-^-U, i = l,2,3, P = PiP2P3 (10) 

leaves Eqs.® invariant, we infer that a point belonging to a solution of Eqs.© lias its counterpart, or twin, at the 
antipodal point and therefore one may visualize the dynamics of L on a sphere with identified antipodal points, that 
is the projective plane. 

Now we are in a position to determine the topological types of the graphs by employing the computer simulation of 
the Eqs.jni) conjunction with the knowledge of the topological types of the fixed points obtained above. It should be 
noted that we must check as to whether the solutions provided by Eqs.JHIl agree with those given by original Eqs.Q, 
see FIG. dlUf) . The phase picture can be obtained by constructing a mesh generated by solutions to Eqs.©, taking 
into account the types of fixed points. The results are illustrated in FIGs(4 — 7). 

We obtain the following topological types of the graphs: 
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FIG. 6; Type I trajectories of the auxiliary system on projective plain; dashed lines are typical trajectories, the solid ones 
separatrixes. Parameters of deformation: ei — 0.02, £2 = 0.03, £3 = 0.04. 




FIG. 7: Type II trajectories of the auxiliary system on the projective plain; dashed lines are typical trajectories, the solid ones 
separatrixes. Parameters of deformation: ei — 0.01, £2 = 0.03, £3 — 0.04. 



Type I FIG01 7 foci and 6 saddles; Ei being subject to the constraints: 

£i£2 - £2e3 +e3ei > 

£i£2 + £263 - e3ei > . (11) 

-6162 + 6263 + SsEi > 

Type II FIGEI 5 foci and 4 saddle points; Ei are not equal to zero, have the same sign, and at least one of eas. Hll|l 
is not true. 

Type III FIGini 3 foci and 2 saddle points; Ei being subject to one of the following constraints: £223 > and £i£2 < 0; 
£3£i > and £2£3 < 0; £i£2 > and £3£i < 0. 

Type IV FIG[7| 2 foci and 1 saddle point; Ei being subject to one of the following constraints: £1 = and £2£3 < 0; 
£2 = and £3£i < 0; £3 = and £i£2 < 0. 



The topological types of the separatrix nets depend on values of the coefficients of the deformation Ei, and generate 
regions I, II, III, IV in the Ei space. Taking into account the homogeneous form of the constraints imposed on Ei, we 
may visualize them on the projective plane corresponding to Ei, see FIG|51 



FIG. 8: Type III trajectories of the auxiliary system on the projective plain; dashed lines are typical trajectories, the solid ones 
separatrixes. Parameters of deformation: ei = —0.02, £2 = 0.03, £3 = 0.04. . 




FIG. 9: Type IV trajectories of the auxiliary system on the projective plain; dashed lines are typical trajectories, the solid ones 
separatrixes. Parameters of deformation: ei = —0.01, £2 = 0.00, £3 = 0.01. 




FIG. 10: Ellipsoid with the main axes: 1.01, 1.02, 1.03; dashed lines are typical trajectories, the solid separatrixes. 
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FIG. 11: Regions of Si corresponding to Types I - IV of the solutions to the auxiliary system, indicated on the projective plain 
represented by a disk with identified antipodal points at the boundary. The white area indicates Type I solutions, the filled and 
the barred ones Type II and Type III. The lines dividing the Type I and Type II regions are subject to Ea. (ll2ll 



It is important that the lines dividing the domains corresponding to types I and II, FIG lllI are given by the 
homogeneous equations 

1. £l£2 - £2^3 + EsEl = (12) 

2. £162 + £2£3 - £3£l = 

3. -£i£2 + £2£3 + £3£l = 

Each equation determines a projective circle on the projective plane of the homogeneous variables £i,£2,£3, see 
FIGsd 

It is important that the solutions to Eqs.(|SJl have the specific feature that the topological type of the graph is 
completely determined by the numbers of foci and saddle points. The dependence of the conformations of the foci 
and the saddles on values of £; is illustrated in FIGEl 

We see that the topological types of the graphs are non-trivial enough, and the graphs realized on the sphere differ 
from those on the projective plain. This circumstance is due to the fact that the graphs on the projective plain are 
obtained from those on the sphere by factoring with transformation l(Tn|) . see FIGlT^ 



III. CONCLUSION 



The key point of the present investigation is the concept of geodetic coil which enables us to cast intuitive geometric 
ideas in analytical form, and relies on constructing an auxiliary hamiltonian system, which can be considered as a 
reduction of the initial problem to that of constructing a graph on projective plane. In analytical terms, one may 
consider it as an asymptotic reduction of the system of equations for geodesies on a deformed sphere to the system 
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FIG. 12: Topological types of the separatrixe nets; Type I — IV. 




FIG. 13: Comparison of the solution to: A. the initial equations for geodesies; B. the averaged equation given by the auxiliary 
system. 



similar to that of the top, with the hamiltonian of the fourth order. The simplification we get in this way, is 
substantial. Indeed, the Hamiltonian system for geodesies could be non-integrable, whereas the auxiliary system is 
totally integrable, and its phase space can be described by a graph that comprises vertices, which correspond to 
stationary solutions, or almost closed geodesies, and edges, which can be visualized as geodetic coils joining them. 
The important thing is that the arguments, of purely analytical and topological nature, which this analysis involves, 
turn out to be helpful for the numerical simulation of the equations for geodetic lines, which admits us to give a 
tangible realization of the visualization problem for geodesies on a surface. 

In fact, this paper is profoundly motivated by the technical means provided by numerical modelling, which has 
allowed us to obtain the final picture of the problem's phase space in terms of the separatrixe graphs. Thus, we feel 
that the approach used in this paper is the symbiosis of the methods of analytical mechanics, computer analysis, and 
topology. The latter is particularly important, giving the conceptual structure for the problem of geodetic lines on the 
deformed sphere. At this point it should be noted that the detailed picture of the phase space reconstruction should 
provide a description of imbedding the invariant torii, which could exist in some regimes, and, perhaps, the domains 
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FIG. 14: Intersection of the conic determining the boundary hne of Type I and the sphere. The projection on the disk gives 
the representation of the hne on the projective plain. 



of chaotic dynamics, peculiar for non-integrable problems, [3. 

The solution obtained in this paper for the geodesies on the deformed sphere does not preclude the chaotic regimes 
of geodetic lines, even though the auxiliary hamiltonian system turns out to be completely integrable. It is worth 
noting that the dynamics of geodesies gives a graphic example of the chaos in riemanniann geometry, 0. This set 
of important problems, lying at the border between nonlinear mechanics, and riemanniann geometry and topology, 
deserves further studying. 
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